Preprint typeset using WT^^ style cmulatcapj v. 10/09/06 



O 
O 

(N 

D 
tin 

in 



> 
in 
m 

O 

o 

:^ 

^: 
6 

■ 



X 



VARIABILITY OF MODERATE LUMINOSITY ACTIVE GALACTIC NUCLEI AT Z=0.36 

Jong-Hak Woo\ Tommaso Treu\ Matthew A. Malkan^, Matthew A. Ferry"', Tony Misch"* 

(Accepted Feb. 5 2007) 

ABSTRACT 

We monitored 13 moderate luminosity active galactic nuclei at z=0.36 to measure flux variability, 
explore feasibility of reverberation mapping, and determine uncertainties on estimating black hole 
mass from single-epoch data. Spectra and images were obtained with approximately weekly cadence 
for up to 4 months, using the KAST spectrograph on the 3-m Shane Telescope. In broad band we 
detect peak-to-peak variations of 9-37% and rms variations of 2-10%. The observed flux variability in 
the g' band (rest-frame 2800-4000A) is consistent with that in the r' band (rest-frame 4000-5200 A), 
but with larger amplitude. However, after correcting for stellar light dilution, using Hubble Space 
Telescope images, we find nuclear variability of 3-24% (rms variation) with similar amplitudes in 
the g' and r' bands within the errors. Intrinsic flux variability of the 11/3 line is also detected at 
the 3-13% level, after accounting for systematic errors on the spectrophotometry. This demonstrates 
that a reverberation mapping campaign beyond the local universe can be carried out with a 3-m 
class telescope, provided that sufficiently long light curves are obtained. Finally, we compare the H/3 
FWHM measured from mean spectra with that measured from single-epoch data, and find no bias 
but an rms scatter of 14%, mostly accounted for by the uncertainty on FWHM measurements. The 
propagated uncertainty on black hole mass estimates, due to the FWHM measurement errors using 
low S/N (10-15 per pixel) single-epoch spectra, is 30%. 
Subject headings: galaxies: active — galaxies: nuclei — quasars: general 



1. INTRODUCTION 

The mass of supermassive black holes (Mbh) is the 
key parameter in understanding the physics of Active 
Galactic Nuclei (AGNs) and the role of black holes in 
galaxy formation and evolution, as indicated in the lo- 
cal universe by the tight relations between Mbh and host 
galaxy properties (Ferrarese & Merritt 2000; Gebhardt et 
al. 2000; Magorrian et al. 1998; Marconi & Hunt 2003; 
Haring & Rix 2004). However, direct dynamical mea- 
surements of Mbh using spatially resolved kinematics of 
stars and gas around the central black hole are limited to 
the local universe, where only about three dozen masses 
have been measured (see Kormendy & Gebhardt 2001; 
Ferrarese & Ford 2005; Tundo et al. 2006). 

Beyond the local universe, Mbh can be estimated for 
active galaxies using the so-called "virial" method. The 
method assumes that the motion of the broad-line re- 
gions (BLR) is dominated by the gravitational potential 
of the central source. Under this assumption, Mbh can 
be determined from the dynamics of the BLR, provided 
that velocity and size are measurable, and that some as- 
sumption on the orbits is made in order to determine the 
virial coefficient. The velocity scale of the BLR can be 
estimated from the width of broad-emission lines, while 
the size (Rblr) can be measured using "reverberation 
mapping" (Blandford & McKee 1982), i.e. determining 
the time lag between continuum changes (presumably 
originating from an accretion flow) and the correspond- 
ing changes in the broad-emission lines. The "virial" 
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assumption has been directly tested for the best studied 
AGN, NGC 5548. Time lags for NGC 5548 have been 
determined using several broad lines of different widths, 
and the inferred distance from the central source shows 
the expected "virial" correlation for a central point mass, 

V cx i?BLR (Peterson & Wandel 2000). 

Since the early attempts of determining the time lag 
between continuum and broad line light curves with 
cross-correlation methods (e.g. Gaskell & Sparke 1986) 
and with more sophisticated response function calcula- 
tions (e.g. Krolik & Done 1995), long-term efforts of re- 
verberation mapping of the local Seyfert galaxies (Wan- 
del et al. 1999) and the low-redshift {z < 0.2) PC quasars 
(Kaspi et al. 2000) have provided clear detections of time 
lags for about three dozens of highly variable AGNs (Pe- 
terson et al. 2004). The measured time lags range from 
a few to ~300 days and the nuclear optical luminosi- 
ties range between 10''^ and 10*^ erg s^^. A Mbh esti- 
mate is then obtained by combining the time-lags with 
measures of the width of the broad lines and a virial 
coefficient. The latter can be obtained for example by 
requiring that AGNs and their host galaxies follow the 
same Mbh -velocity dispersion (cr; or bulge luminosity) 
relation as quiescent galaxies (Onken et al. 2004; Greene 
& Ho 2006, Labita et al. 2006). Mbh estimates via re- 
verberation mapping are believed to be accurate within 
a factor of ~3 (Bentz et al. 2006b; Onken et al. 2004; 
however, see Krolik 2001). 

Reverberation mapping based masses have also pro- 
vided a widely used empirical method to estimate Mbh 
in distant AGNs from single-epoch data, bypassing the 
need for observationally expensive time series. This em- 
pirically calibrated photo-ionization method (Wandel et 
al. 1999) is based on the observed correlation between 
the BLR size (time lag) and optical luminosity at 5100A. 
The slope of the power law relation Rblr, oc ^5100 with 
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a ~ 0.5 is similar to that expected by photoionization 
models of the BLR (Wandel et al. 1999; Kaspi et al. 
2000; Bentz et al. 2006a). This empirical relation has 
been widely used to determine Mbh for large samples 
of broad-line AGNs, where long-term monitoring efforts 
would have been prohibitive (e.g. Woo & Urry 2002; 
McLure et al. 2004). The uncertainty of Mbh estimates 
based on this empirical relation is believed to be a factor 
of 3-4 (Vestergaard 2002, 2006), and it is due to sev- 
eral factors: i) intrinsic uncertainties in the Mbh of lo- 
cal calibrators; ii) intrinsic scatter of the size- luminosity 
relation; iii) differences between single-epoch measure- 
ments of optical luminosity and width, typically available 
for high-redshift studies, and those obtained from mean 
or rms spectra obtained from multiple epochs, generally 
used in the reverberation mapping analysis. 

Determining Mbh at high redshift is crucial to make 
progress in a number of outstanding scientific issues, e.g., 
black hole demographics (Yu & Tremaine 2002), accre- 
tion mechanisms (Koratkar & Blaes 1999), connections 
between galaxy formation and evolution and AGN feed- 
back (Croton 2006; Hopkins et al. 2006), origin of the 
Mbh-c and other scaling relations (e.g. Silk & Rees 
1998; Kauffmann & Haehnelt 2000). Studies of AGN 
variability at high-redshift and possibly direct reverber- 
ation mapping determinations of time lags would be ex- 
tremely valuable to validate the photoionization method 
and understand its uncertainty. However, reverberation 
mapping is challenging for bright quasars at high-redshift 
because it requires a longer time base line due to the time 
dilation (l + z) effect. Moreover, high luminosity objects 
(L5100 > 10'*^ erg s~^) seem to have smaller amplitudes of 
continuum variability compared to Seyfert galaxies (vari- 
ability amplitude-luminosity anticorrelation; Cristiani et 
al. 1997; Vanden Berk et al. 2004). For example, af- 
ter six year of monitoring, Kaspi et al. (2006) detected 
no variability of the Lya line for 6 high-redshift (z=2.2- 
3.2) quasars while they detected C IV A1550 variability 
(see also Wilhite et al. 2006). These high luminosity 
quasars show a factor of two smaller continuum variabil- 
ity than that of the lower luminosity PG quasars moni- 
tored over a comparable rest-frame period. Lower lumi- 
nosity type-1 AGNs (Seyfert Is), although fainter, pro- 
vide several advantages from the point of view of a mon- 
itoring campaign. First and foremost the expected time 
lags are shorter (order of weeks to a few months) , raising 
the hopes of detecting a lag in only one season, with- 
out having to worry about objects observability. Second, 
variability is expected to be more pronounced since the 
amplitude of variability seems to be anticorrelated with 
luminosity. 

In this paper we report on our time-domain study of 
13 moderate luminosity (L5100 ~ lO'*'* erg s~^) AGNs at 
z ^0.36. The campaign was carried out at the Lick 3-m 
telescope from May 2004 to November 2004. The sample 
is drawn from our sample of Seyferts, for which stellar ve- 
locity dispersion (Treu, Malkan & Blandford 2004; Woo 
et al. 2006; hereafter TMB04, W06) and bulge luminosity 
and radius (Treu et al. 2006 in preparation) have been 
determined from Keck spectroscopy and Hubble Space 
Telescope imaging in order to investigate the cosmic evo- 
lution of the Mbh-c and Mnn-bulge luminosity relations. 
The goals of this paper are to: i) study the variability of 
distant seyferts to explore the feasibility of reverberation 



mapping in the distant universe; ii) study the variability 
of the H/3 line width to determine the contribution to the 
systematic uncertainty of the Mbh estimated from single 
epoch spectra as opposed to rms or mean spectra. The 
paper is organized as follows. In § 2 we describe observa- 
tions and data reduction. In § 3 we describe continuum 
and line flux variability. In § 4 we compare line width 
measurements from single-epoch data. In § 5 we discuss 
our results and their implications. A standard cosmol- 
ogy is assumed where necessary (Ho=70 kms^^Mpc^^, 
= 0.3 and Qa = 0.7). 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Experiment Design and Sample Selection 

In order to increase our chances to detect a time lag 
in a single season, we focused our efforts on moderate 
luminosity AGNs, for which the time lag is predicted to 
be of order a few weeks in the rest frame, based on the 
size-luminosity relation (Kaspi et al. 2005). A contin- 
uos coverage is critical for a monitoring campaign, but 
the Lick observatory does not offer service mode obser- 
vations. Thus, we decided to carry out our campaign 
over the summer months, when the typical weather pat- 
terns at Mt. Hamilton indicate that a relatively large 
fraction of time is useful for observations (above 80% 
between June and September^). Note that this program 
does not require photometric conditions, nor particularly 
good seeing. Of course, the drawback of observing in the 
summer is the short duration of the nights. Given the 
expected time lags, we planned observations with weekly 
cadence (corresponding approximately to 5 days in the 
rest frame). The Mt. Hamilton and Lick Observatory 
staff helped enormously by scheduling with regular ca- 
dence, except for the obviously longer gaps during full 
moon periods when the background was prohibitive for 
optical observations. 

The targets are drawn from a parent sample of 30 
broad-line AGNs at z=0.365±0.010, initially selected 
from the Sloan Digital Sky Survey (SDSS) Data Release 
2 (Abazajian et al. 2004) to investigate the relation be- 
tween Mbh and their host galaxy properties (for details 
of sample selection, see W06). A subsample of 13 targets 
was selected based primarily on observability during the 
summer months, and secondarily on the availability of 
ancillary Keck and Hubble data (Table 1). 

2.2. Observations 

All observations were performed with the Shane 3-m 
telescope at the Lick Observatory between May 10 and 
November 19. Table 2 shows the journal of observations. 
The weather was generally bad in May and after the 
end of September, effectively cutting our campaign to 
4 months. A fraction of time was also lost to techni- 
cal problems. When we could open the dome, observing 
conditions ranged from photometric to thick cirrus, with 
a typical seeing ~2" . In total, more than 100 individual 
spectra and images were taken on 20 nights, allowing us 
to measure variability of continuum and line fluxes. 

We used the Kast double spectrograph with the 600 
line mm^-'^ grism centered at ~ 4300 A to obtain blue 
(3300-5500 A) spectra, and with the 600 line mm^^ grat- 

^ See |http://mthamiIton.ucoIick.org/techdocs/MH_weather/obstats/prcnt-hrs 
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ing centered at ~ 7500A to obtain red spectra (6300- 
9IOOA). The pixel scale corresponds to 1.83 AxO'.'S in 
the blue and 2.32AxO'.'8 in the red. The blue spectra 
cover the Mgll line (rest wavelength ~ 2798A) while the 
red spectra cover the H/3 line (rest wavelength 4861 A) 
and extend in most cases to the red to include the Ha 
line. With the 2" wide slit , the spectral resolution (gaus- 
sian dispersion) measured from arc lines was ~ 180 km 
s^^ around the Mgll line and ~ 100 km s~^ around the 
H/3 line. Typical exposure time for each object was 2700- 
3000 second or 2 x 1800 second, yielding S/N - 10 per 
pixel on the continuum of single-epoch spectra. 

Internal flat fields for the red spectra were obtained 
at each target position to correct pixel-to-pixel variation 
and the fringing pattern of the red CCD. For the blue 
spectra, internal flats were taken in the afternoon. A set 
of AOV stars - selected from the Hipparcos catalog to 
be close in the sky to our targets - was used to correct 
the A and B-band atmospheric absorption features and 
to perform secondary flux calibration (see TMB04 and 
W06 for details) . Spectrophotometric standards were ob- 
served for flux calibrations. 

For spectroscopic target acquisition and to obtain 
broad band imaging time-series, we used the direct imag- 
ing capability of the Kast spectrograph. Before every 
spectroscopic exposure, we obtained images in the g' and 
r' bands over the unvignetted field of view of approxi- 
mately 2 arcminute diameter (0.8 arcsec pixel"^). The 
field of view is large enough that for most objects at least 
a couple of stars brighter than the target AGN were avail- 
able for differential photometry. Typical exposure time 
for imaging was 120 second on the red side with the r' 
filter and 150 second on the blue side with the g' filter. 

2.3. Data Reduction 

We performed the standard data reduction includ- 
ing flat fielding, wavelength calibration, spectral extrac- 
tion, and flux calibration using a series of IRAF scripts. 
Cosmic rays were removed from each exposure using 
the Laplacian cosmic-ray identification software (van 
Dokkum 2001). Sky emission lines were used for wave- 
length calibration in the red, supplemented by arc lamps 
in the blue. One-dimensional spectra were extracted for 
maximal S/N with a typical extraction radius of 3 pix- 
els, corresponding to ^ 4.8". Thirteen low quality spec- 
tra with low S/N (mostly due to bad weather) or severe 
fringing effects (due to imperfect flat-fielding) were re- 
moved from the high quality sample to be analyzed for 
variability. The rest-frame mean and rms spectra for all 
13 AGNs are shown in Figure 1, together with the aver- 
age noise level. Broad Mgll, H/3, Ha and narrow OIII 
lines are clearly visible in all the average spectra. In most 
cases, the H/3 line is not visible in the rms spectra since 
the variation over the observed epochs is smaller than 
the average uncertainty per pixel. The Ha line is clearly 
visible in some of the rms spectra. 

3. FLUX VARIABILITY 

Variability is one of the main characteristics of AGNs, 
and is the property that is relied upon for reverberation 
mapping studies. However, the amplitude and the spec- 
tral shape change is not well studied for general Seyfert 
galaxies and quasars. In this section, we first present 



optical continuum variability (§ I3.ip as determined by 
broad band g' and r' photometry, using available HST 
photometry to separate the variable nuclear component 
from the constant stellar component (§ 13. 2p . Then, in 
§ 13.31 we study the H/3 line fiux variability from the spec- 
troscopic analysis. We focus on the H/3 line since this is 
the most commonly used line in this redshift range, and 
the proximity of the narrow OIII lines provides a robust 
relative spectrophotometric calibration. 

3.1. Continuum flux variability from Lick photometry 

In order to measure continuum variability, we per- 
formed differential photometry using stars in the field 
of each AGN with both g' and r' band images. Three 
to five stars around the target AGN were typically used 
to estimate relative flux variations of the AGN. We ex- 
cluded stars at the edge of the images, where vignetting 
is problematic. Aperture photometry of individual stars 
and AGNs was performed using the phot task in IRAF, 
with a 3 pixel (2.4") aperture radius for flux measure- 
ments and with an annulus between 8 and 12 pixels for 
sky subtraction. The optimal size of the aperture and sky 
subtraction annulus were determined based on extensive 
tests with images of various quality seeing. Typical er- 
rors in the aperture photometry are 1-3% for the 17-18th 
magnitude reference stars. 

Magnitude differences between AGN and field stars 
were calculated and then normalized to zero by subtract- 
ing the mean difference over all observed epochs. Pho- 
tometric errors of the target AGN and each star were 
added in quadrature, yielding total errors on differential 
photometry in the range 2-4%. Normalized magnitude 
differences obtained with different stars are consistent 
within the errors, providing a good sanity check on the 
analysis. The final light curve and errors were produced 
by averaging normalized light curves and errors for indi- 
vidual reference stars. 

We measured the continuum fiux variability for 7 ob- 
jects, namely S04, S05, S06, S24, S27, S40, and S99, with 
more than 5 epochs of reliable observations (see Fig. 2, 
3 and 4). A continuum broad band light curve for S08 
could not be obtained due to the lack of suitable reference 
stars in the field. The observed flux variability in the 
g' band, corresponding to rest-frame wavelength 2800- 
4OOOA, is consistent with that of the r' band (rest-frame 
wavelength 4000-5200A), but with larger amplitude. In 
the g' band, peak-to-peak variation is 12-37% and rms 
variation is 2-10% while in the r' band peak-to-peak vari- 
ation is 9-32% and rms variation is 2-9% (see Table 3). 
This amount of variability on week to month timescales is 
consistent with previous studies of local Seyfert 1 galax- 
ies and quasars (e.g. Webb & Malkan 2000; Klimek et 
al. 2004) and high-redshift quasars (Kaspi et al. 2003). 

3.2. Continuum nuclear flux and color variability 

Stellar light contributions from host galaxies to the to- 
tal flux is often signiflcant, especially for Seyfert galaxies 
(Malkan & Filippenko 1983; Bentz et al. 2006a). This 
leads to reduced variability in the integrated spectrum 
and hampers studying color variability of the nucleus 
(Hawkins 2003). 

To investigate nuclear flux and color variability, we 
need to correct for the host galaxy contamination. This 
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can be achieved using the AGN-galaxy decomposition 
analysis obtained from HST-ACS F775W (i') band im- 
ages (GO-10216; PI Treu; Treu et al. 2006, in prepara- 
tion), and adopting our single-epoch measurement of the 
nuclear light fraction as our best estimate of the aver- 
age nuclear light fraction. First, we derived the stellar 
fraction for the g' and r' bands using: 



(1) 



Here, fm is the stellar fraction in each band, li is the 
stellar fraction in the i' band determined from the HST 
imaging analysis, (m-i)jotai is the observed color from 
the SDSS photometry (DR4), and (m-i)^^; is the stellar 
light color, g-r = 1.65 ± 0.15 and r - i = 0.60 ± 0.06, 
as estimated from population synthesis models (Bruzual 
& Chariot 2003). The uncertainty on the stellar colors 
represents the spread obtained from models with ages 
ranging from 1 to 7 Gyr and metallicities from 0.4 to 
2.5 solar. Then, we derived intrinsic flux variability by 
correcting the host galaxy contamination with: 



Rc — 



R-f 



(2) 



Here, f is the stellar fraction in each band, R is the ob- 
served variability (i? = Ftotai/ < Ftotai >), and Rc is 
the corrected variability {Rc = Fagn / < Fagn >)• The 
mean correction on the variability of 6 objects with HST 
images is 0.016 magnitude in the g' band and 0.057 mag- 
nitude in the r' band. After correcting host galaxy con- 
tamination, we find 3-14% rms variation in the g' band 
and 5-24% in the r' band, somewhat larger in the red for 
a few objects (see Table 3). 

Figure 3 compares the continuum flux variability in the 
g' and r' bands with and without host galaxy correction. 
Uncorrected flux variability shows higher amplitude in 
the g' band than in the r' band {Ag'/Ar' = IMt^Ts)' 
as expected because stellar light dilution is more promi- 
nent in the red. However the variations become much 
more similar when the host galaxy contaminations are 
corrected for {Ag'/Ar' = 0.71 ± o'mt"-]/^^ where errors 
include a random component and a systematic error from 
uncertainties in stellar colors). This suggests that the 
spectral shape in the rest-frame 2800-5200A does not 
significantly change on weekly time scales, similar to the 
findings of Wilhite et al. (2005) that shows an ensem- 
ble quasar spectrum is bluer in bright phases only at 
rest wavelength < 2500A. The slightly redder colors in 
bright phases are probably caused by the lack of (or much 
lower) variability of emission lines, which are included in 
the wavelength coverage of the brodbands (MgH in blue, 
and H/3 and [O III] in red), because the relative emission 
line contribution to the continuum flux is lower in bright 
phases, producing redder photometric colors (Wilhite et 
al. 2005). 

3.3. Flux variability of the Hj3 line 

Emission line flux variability is more difficult to mea- 
sure than continuum variability, since the uncertainty in 
the flux calibration of the ground-based spectroscopy is 
typically larger than 10%, particularly for faint sources. 
In the case of AGNs, however, narrow emission lines such 
as [O III] can be used as internal flux calibrators, since 
these lines originate from extended low-density regions 



- probably more than 100 times further from the accre- 
tion disk than the BLR (Bennert et al. 2006) and 
do not vary on short time scales as broad-emission lines 
do. Thus, flux variations of broad lines can be measured 
by normalizing each spectrum to the constant narrow 
line fluxes, eliminating uncertainties related to slit loss 
effects, sky transparency, airmass correction, etc. An 
additional advantage of the specific case of the doublet 
[O III]AA4959,5007 is that the flux ratio between the two 
components does not vary (Bachall et al. 2004, Dimitri- 
jevic et al. 2007) and therefore the rms of the measured 
ratio provides a robust estimate of the residual uncer- 
tainties on the normalized flux. 

For 8 objects (S04, S05, S06, SOS, S24, S27, S40, S99) 
with more than 5 reliable single-epoch spectra, we mea- 
sure H/3 line flux variations using [O III] A 5007 line as 
an internal flux calibrator as described below (the pro- 
cedure is similar to that adopted by local reverberation 
mapping studies, e.g. Peterson et al. 2002). We also 
measure the [O III] A 4959 to [O III] A 5007 line ratio to 
help optimize the wavelength windows used for line flux 
determination^ . 

First, we subtract the continuum under H/? and 
[O III] lines using the average flux around 4700A and 
5100A with ~ 50A windows. We linearly interpolate 
these two continuum points to define the underlying con- 
tinuum. Second, we measure the H/3 and [O III] doublet 
line fluxes by integrating the flux over each line in the op- 
timal windows. Third, the H/3 and [O III] A 4959 fluxes 
are normalized by the [O III] A 5007 flux. 

Before we can proceed to determine the intrinsic H/3 
line flux variability, we need to account for residual sys- 
tematic errors in the line ratios (csyg). We use the fact 
that the flux ratio of [O III] A 5007 and 4959 is constant 
to obtain a conservative estimate of agys- This is done 
by requiring that a constant ratio be an acceptable fit to 
the data as measured by the statistic: 



E 



[R{i)- < R{i) >Y 



N -I. 



(3) 



Here, R(i) is the flux ratio between [O III] A4959 and 
A5007 at each observed epoch, a{i) is the error in the 
flux ratio, (Jsys is the systematic error, and N is the num- 
ber of observed epochs. The inferred systematic errors 
are in the range 2-15%. We note that this is likely to 
be an upper limit to the systematic errors as A4959 is 
weaker than H/3, and narrow lines fluxes are more sen- 
sitive to uncertainties on wavelength calibration and to 
poor wavelength sampling than broad lines. 

Taking random and systematic errors into account, we 
can now estimate the H/'? intrinsic line variability (Ta- 
ble 3). We flnd rms variability in the range 3-13%. We 
caution that this should be considered a lower limit due 

® The optimal wavelength windows for narrow line flux determi- 
nation are a tradeofi' between the advantages of a small window — 
high signal-to-noise per pixel, little effects of continuum subtrac- 
tion uncertainties — and the disadvantages — sensitivity to uncer- 
tainties in wavelength calibration as well as to pixelization in the 
wavelength domain. To find the optimal window we repeated the 
procedure described below for a range of wavelength regions and 
picked the one that minimizes the scatter in the [O III] doublet 
line flux ratio. The choice of the wavelength region for the H/3 line 
is far less critical since the decline in S/N is more gradual and the 
spectrograph's pixels are small compared to the line width. 
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to the fact that our estimate of systematic errors is con- 
servative. For 3 objects we did not detect significant 
variabihty. 

3.4. Prospects for reverberation mapping 

Wc compare flux variations between continuum and 
the H/3 line flux in Figure 4. The detected H/3 line vari- 
ability is smaller than the rms variations of the contin- 
uum flux (a factor of 2 smaller on average after host 
galaxy correction), similar to other works (e.g. Rosen- 
blatt & Malkan 1992). Errors on the H/3 flux measure- 
ments are larger than those on the continuum flux mea- 
surements because the S/N of the spectra is significantly 
lower than that of broad band images. The light curve 
of the H/? flux is qualitatively different from that of con- 
tinuum flux, possibly indic;ating a time lag. However, 
wc do not detect a clciar lag between continuum and H/3 
light curves. This could be due to the fact that the time 
base line of our campaign is relatively short (4 months 
at maximum) although the observed variation amplitude 
(5-10% on H/3) is comparable to that of the time-lag mea- 
sured samples in the literature. To illustrate this, in Fig- 
ure 5 we compared the monitoring time base lines in the 
rest-frame of each object with the variation amplitude 
for time-lag detected/undetected samples from the liter- 
ature. In general a few hundred days of time base lines 
were required to detect a reliable lag. For high luminos- 
ity PG quasars, several years of base lines were used to 
get a time-lag. Out of 3 objects in our sample (open 
circles) with the longest time base line (~100 days in 
the rest-frame), one object does not show H/3 variability 
while the other two objects show comparable variability 
as the Scyfcrt galaxies with a time lag detection (flUcd 
squares). In contrast, rcst-framc time base line for our 
sample is at best in the borderline of the regime, where 
a reliable time lag can be measured, suggesting that we 
may detect a time lag if a factor of ^2 longer time base 
line can be provided. 

4. H/3 LINE WIDTH FROM SINGLE-EPOCH DATA 

The BLR velocity scale is one of the two ingredients in 
estimating virial Mbh- In reverberation mapping Mbh 
estimates, the FWHM of Balmer lines is generally mea- 
sured from mean or rms spectra averaged over many ob- 
served epochs (Kaspi et al. 2000). In contrast, for most 
of Mbh studies using the size-luminosity relation, rela- 
tively low S/N data from a single-epoch observation, e.g. 
SDSS quasar spectra, are used. This could introduce ad- 
ditional uncertainty to virial Mbh estimates, if the width 
of broad lines varies significantly or the random error on 
FWHM measurements on single-epoch data is significant. 
Here, we investigate the uncertainty and variation of the 
H/3 line width by comparing FWHM from a mean spec- 
trum with that from each spectrum of various observed 
epochs. 

The FWHM of the H/3 line is measured in several steps. 

First, we subtract continuum under the H/3 line by iden- 
tifying the continuum levels of each side of the H/3 line 
and interpolating between them. Second, we subtract 
[O III] A4959 by dividing [O III] A5007 by 3 and blueshift- 
ing. Third, we subtract the narrow component of H/3 by 
rescaling [O III] A5007 and blueshifting it. Since the 
FWHM measurements are sensitive to the ratio of nar- 
row H/3 to [O III] (up to ~10%), we used the fixed ratio 



determined from the high S/N Keck spectra (sec Table 
3). Fourth, wc fit the broad component of H/3 with the 
Gaussian-Hermite polynomials up to 5-6 orders, depend- 
ing on the asymmetry of the H/3 line profile, and measure 
the FWHM of the model (Figure 6). This is necessary 
because our single-epoch spectra have a typical S/N of 
10-15, too low to precisely define wavelengths at the half- 
maximum. 

Figure 7 compares the H/3 FWHM measured from 
mean spectra with that measured from each single-epoch 
data for 9 objects with more than 3 reliable spectra. 899 

is excluded in this analysis since the H/3 line has such a 
weak and very broad profile that we could not get mean- 
ingful measurements from single-epoch data. With 59 
single-epoch data with S/N > 10, we find no system- 
atic bias with an average ratio of 1.03 ± 0.02 between 
single-epoch FWHM and FWHM from mean spectra. 
The rms scatter around FWHM from mean spectra is 
^ 14%, similar to 15-20% scatter found by Vestcrgaard 
(2002), who compared one single-epoch FWHM with 
FWHM from mean spectra for each of 18 AGNs. This 
scatter is a combination of actual width variation over 
monthly timescales, and random and systematic errors 
on single-epoch measurements. For example, Rosenblatt 
et al. (1992) measured H/3 FWHM using 12 single-epoch 
spectra for 12 local Seyfert galaxies and found 7-26% rms 
scatter (mean rms scatter 15%), which were comparable 
to their FWHM measurement uncertainty 

To understand intrinsic width variation and measure- 
ment errors, we measured H/3 FWHM for 8 objects (filled 
squares), for which high S/N (50-100) Keck single-epoch 
data are available (W06). The random errors on the 
FWHM measurements from the Keck spectra are virtu- 
ally negligible. The rms scatter between the H/3 FWHM 
from Lick mean spectra and that from Keck single-epoch 
spectra is 7% (Figure 7). We note that this width vari- 
ation is much smaller than continuum flux variation as 
expected from the size-luminosity relation in which log 
FWHM scales with 1/4 x log L for given Mbh- If we take 
this scatter as a measure of the intrinsic width variabil- 
ity, then we can consider errors on the Lick single-epoch 
FWHM to be 12% (0.05 dex), by subtracting 7% (0.03 
dex) variation from 14% (0.057 dex) scatter in quadra- 
ture. Thus, most of the variation in the single-epoch 
FWHM measurements seems to be caused by the mea- 
surement uncertainty. Consequently, we do not find any 
correlation between variations in FWHM and continuum 
or line flux since any intrinsic trend is swamped by mca- 
siircmcnt errors. If we takc> 14% scattcir as the uncer- 
tainty on FWHM measurements, the propagated error in 
Mbh estimates based on single-epoch data is 30%. This 
estimate should be applied with caution to other single- 
epoch data sets since the 30% propagated error may not 
be relevant because of the difference in S/N, the method 
of FWHM determination, variability, and other factors. 

5. DISCUSSION AND CONCLUSIONS 

We detected variability of moderate luminosity 
(L5100 ^ 10'"' erg s-i) AGNs at z=0.36. Continuum flux 
variability in the rest-frame wavelength 2800-4000A and 
4OOO-52OOA is detected at the 2-10% level (rms vari- 
ation). These measurements are consistent with other 
variability studies of local Seyfert galaxies and higher 
redshift quasars (e.g. Webb & Malkan 2000; Kaspi et al. 



2000). After correcting for stellar light contamination us- 
ing HST-ACS images, we find that the intrinsic nuclear 
continuum variability is 3-25%. Blue and red continuum 
variability correlates, indicating that the spectral shape 
does not significantly change on weekly time scales. 

We also detected intrinsic H/3 line flux variability at the 
level of 3-13%, on average a factor of 2 lower than the 
continuum flux variability. Light curves of the contin- 
uum flux and the H/3 line flux are qualitatively different, 
perhaps due to a time lag. However, the time base line of 
our campaign seems to be too short to detect a reliable 
lag. We may detect a time lag if a longer time base line 
can be provided. With modern scheduling techniques 
and perhaps a robotic telescope, light curves could be 
extended to 8 months per year. 

Finally we study variability of the H/3 line width to 
estimate uncertainties in Mbh estimates from single- 
epoch data. We compare the H/3 FWHM measured from 
high S/N mean spectrum with that measured from each 
single-epoch spectrum. We flnd 14% scatter around the 
one-to-one relationship, most of which comes from ran- 
dom errors on FWHM measurements using single-epoch 
spectra. Mbh estimated from single-epoch spectra with 
typical S/N of 10-15 has 30% uncertainty due to the line 
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TABLE 1 
Targets 



— 

(1) Name 




(2) z 




(3) RA (J2000) 


(4) DEC (J2000) 




(5) r' 


804 n2in2n 50-064645 0) 


3580 


21 02 11.51 


-06 46 45 03 


18.75 


S05 (J210451. 83-071209.4) 


0.3531 


21 04 51.85 


-07 12 09.45 


18.43 


S06 (J212034.18-064122.2) 


0.3689 


21 20 34.19 


-06 41 22.24 


18.53 


SOS {J235953.44-093655.6) 


0.3591 


23 59 53.44 


-09 36 55.53 


18.61 


S09 {J005916.10+153816.0) 


0.3548 


00 59 16.11 


+ 15 38 16.08 


18.33 


S12 {J021340.59-f 134756.0) 


0.3575 


02 13 40.60 


+ 13 47 56.06 


18.18 


S23 {J140016.65-010822.1) 


0.3515 


14 00 16.66 


-01 08 22.19 


18.22 


S24 (J140034.70+004733.3) 


0.3621 


14 00 34.71 


+00 47 33.48 


18.43 


S26 (J152922. 24+5928.54.5) 


0.3691 


15 29 22.26 


+59 28 54.56 


18.93 


S27 (J153651. 28+541442. 6) 


0.3667 


15 36 51.28 


+54 14 42.71 


18.87 


S28 {J161156.29+451611.0) 


0.3682 


16 11 56.30 


+45 16 11.04 


18.78 


S40 (J012655.82+153357.8) 


0.3749 


1 26 55.82 


+ 15 33 57.87 


18.62 


S99 (MS 1558.3+4138) 


0.3690 


16 GO 02.80 


+41 30 27.00 


18.78 



Note. — Col. (1): Target ID. Col. (2): Redshift. Col. (3): RA. Col. (4): DEC. Col. (5): Extinction corrected r' AB magnitude from 
SDSS photometry. 



TABLE 2 
Journal of observations 



Run Date 


Targets 


Seeing 
aJTcsec 


Conditions 


(1) (2) 


(3) 


(4) 


(5) 


1 2004 May 10 


S26,S99 


2.4-3 


clear / fog,humidity 


2 2004 May 18 


S05,S23,S24,S27,S99 


1.7-1.8 


clear 


3 2004 May 25 


S23,S24,S26,S27,S28,S99 


2.0-2.2 


photometric 


4 2004 Jun 1 


S05,S24,S26,S27,S28,S99 


1.8-2.2 


photometric 


5 2004 Jun 11 


S05,S24,S26,S27,S28,S99 


1.8-2.4 


cirrus 


6 2004 Jun 20 


S05,S24,S26,S27,S28,S99 


2.1-2.8 


cirrus 


7 2004 Jun 29 


S05,S24,S26,S27,S28,S99 


2.2-2.9 


clear 


8 2004 Jul 9 


824,827,899 


1.5-1.6 


clear/instrument failure 


9 2004 Jul 14 


805,S24,S27,S28,S99 


2.2-2.4 


clear 


10 2004 Jul 18 


S04,S05,S06,S27 


2-3 


cirrus 


11 2004 Aug 6 


S04,S05,S06,S08,S09,S27,S99 


1.7-2. 


photometric/CCD problem 


12 2004 Aug 14 


S04,S05,S06,S08,S09,S27,S99 


2.2-3 


clear /cirrus 


13 2004 Aug 24 


S04,S05,S06,S08,S09,S27,S99 


1.7-2 


cirrus 


14 2004 Sep 8 


S04,S05,S06,S08,S09,S12,S27,S40.S99 


1.3-1.7 


clear /guider problem 


15 2004 Sep 15 


804,805, S06,S08,S09,S12,S27,S40,S99 


1.4-1.8 


photometric 


16 2004 Sep 22 


804,S05,S08, 309,812, S27,S40,S99 


1.3-2.4 


cirrus 


17 2004 Sep 30 


804,805,806,808,809,827,840 


1.4-1.8 


cirrus 


18 2004 Oct 11 


804,805,806,808,809,812,827,840 


1.8-2.2 


photometric/instrument problem 


19 2004 Oct 20 


809,840 


2-3 


cloud, humidity 


20 2004 Nov 3 






rain 


21 2004 Nov 19 


812,840 


2.5-2.7 


clear 



Note. — Col. (1); Observing run. Col. (2); Observing date. Col. (3): targets. Col. (4): Seeing (FWHM) measured from r' band 
direct imaging. Col. (5): Conditions/Notes. 
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TABLE 3 

Variability 



-L til t 


of cpoclis 


TTICIXq' 


TTflS q' 








TTflS J.I 




Tins J.' f. 




O'sys. 




II ft /[n TTTl 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


S04 


8/8/9 


0.338 


0.098 


0.026 


0.118 


0.096 


0.022 


0.024 


0.058 


0.102 


0.068 


0.060 


0.07 


S05 


14/13/10 


0.121 


0.026 


0.026 


0.029 


0.111 


0.028 


0.024 


0.049 


0.084 


0.051 


0.081 


0.10 


S06 


9/8 /5 


0.217 


0.061 


0.064 


0.079 


0.208 


0.042 


0.043 


0.127 


0.014 






0.10 


SOS 


0/0/8 


















0.200 


0.162 




0.10 


S24 


6/7/6 


0.245 


0.078 


0.044 


0.098 


0.300 


0.099 


0.037 


0.243 


0.111 


0.053 


0.049 


0.08 


S27 


16/17/15 


0.289 


0.078 


0.045 


0.090 


0.210 


0.034 


0.033 


0.075 


0.165 


0.107 




0.18 


S40i 


5/5/5 


0.147 


0.032 


0.057 




0.100 




0.054 




0.035 


0.016 


0.029 


0.07 


S99 


14/15/11 


0.338 


0.072 


0.050 


0.136 


0.197 


0.040 


0.028 


0.193 


0.066 


0.020 


0.130 


0.10 



References. — 1) no HST image available. 

Note. — Col. (1): Target ID. Col. (2): Number of flux points in the g' and r' band images and H/3 line flux. Col. (3): peaJc to peak 
variation in magnitude in g' Col. (4): rms variation in g' Col. (5); mean error in g' Col. (6): rms variation in g' after stellar contribution 
correction. Col. (7): peak to peak variation in magnitude in r' Col. (8): rms variation in r' Col. (9): mean error in r' Col. (10): rms 
variation in r' after stellar contribution correction. Col. (11): rms scatter of the [O III] A5007 to [O III] A4959 flux ratio. Col. (12): 
Systematic errors on estimated as described in Section 3.3 Col. (13) Intrinsic rms H/3 flux variation in magnitudes after correcting 
systematic errors. Col. (14) narrow H/3 component to [O III] A5007 flux ratio, used to remove narrow H/3 component before measuring 
broad H/3 FWHM. 
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Fig. 1. — Mean (black) and rms (blue) spectra of all observed AGNs, covering Mgll, H/9, and Ho lines, with the average noise level 
(red). The number of combined spectra is designated for each objects. For S06 and S99, blue spectra are excluded since Mgll line is not 
well defined due to the much lower S/N ratios. 
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Fig. 1. — Continued 
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Fig. 2. — Examples of light curves in the g' and r' bands. Differential photometry was obtained by comparison with nearby stars in the 
imaging field as described in Section 3.1. For clarity, g' band light curves (blue circles) are offset by 0.2 magnitudes. 




Fig. 3. — Comparison of variability in g' and r' bands. Left: observed flux variability. Individual AGNs are denoted with different colors 
(green: S04, magenta: 305, yellow: S06, black: S24, blue: S27, red: S99). Note that variations in g' and r' band are correlated, with a 
larger amplitude in the blue band. The solid line is the best fit slope, 1.84^Q'jg, significantly larger than unity (dashed line), suggesting 
that spectra look bluer when brighter. Right: Nuclear flux variability after correcting for the stellar contribution using HST-ACS images. 
Note that the amplitude of the nuclear variations is larger and the plotting range is larger than in the left panel. The best flt slope is now 
0.71 ± Omto.ll, where errors include a random component and a systematic error from uncertainties in stellar colors. This is marginally 

consistent with unity, suggesting approximately achromatic intrinsic luminosity variations over this wavelength range 2800-5200A. Typical 
errors are plotted at the bottom-right of each panel. 
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Fig. 4. — Light curves of the continuum flux and the H/3 Une flux for SOS and S27. For each object the continuum flux variability is 
shown in the g' band (top) and the r' band (middle). The H/3 line flux light curve is shown in the bottom panel. Note that the amplitude 
of the variations is similar to that of the broad band. 
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Fig. 5. — H/3 variability vs. rest-frame time base line for various samples. Time lag detected Seyfcrts (filled squares) and PG quasars 
(filled triangles) are from the compilation of Peterson et al. (2004). AGNs with no or unreliable time lag detection are denoted with open 
symbols (circles: our Seyferts at z=0.4; triangles: 5 AGNs from Lovers of AGNs campaign, Jackson et al. 1992; Dietrich et al. 1994; Stirpe 
et al. 1994; Erkens et al. 1995). Time base line for time lag measured sample is generally more than a few hundreds days in the rest-frame. 
In contrast, time base line for our sample is relatively short (<100 days) although variability of the H/3 line flux is comparable to that of 
the time-lag measured Seyferts. 
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Rest wavelength (A) 

Fig. 6. — Example of the H/3 FWHM measurement using Gaussian-Hermite polynomials. After removing the narrow H/3 and [O III] 
A4959 with a rescaled and blueshifted [O III A5007 line profile (green), Gaussian-Hermitian models (red) are used to fit the broad H/3 line 
(blue). The FWHM of the line is then measured on the model fit. 
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Fig. 7. — Comparison of the H/3 FWHM measured from mean spectra with those measured from various single-epoch spectra. Single- 
epoch FWHMs are consistent with FWHM from mean spectra with an average ratio of 1.03 ± 0.02. The rms scatter around FWHM 
from mean spectra is 14%, indicating uncertainties of Mbh estimates based on single-epoch data can be 30% due to the random errors in 
measuring FWHM. For comparison, we show FWHM measurements using high S/N Keck data (filled squares), which are consistent with 
FWHM of mean spectra with a rms scatter ~0.03 dex, indicating intrinsic FWHM variation is ~7%. 



